function dYdt = nonlin_state(t,Y,A)

%Cross terms

ind = 1:length(Y);
shift_ind = ind([2:end 1]);

Ycross = Y.*Y(shift_ind);

%Linear transformation of nonlinear states

dYdt = A*Ycross;